Nuclear model effects in Charged Current neutrino-nucleus quasielastic scattering 
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The quasielastic scattering of muon neutrinos on oxygen 16 is studied for neutrino energies between 
200 MeV and 1 GeV using a relativistic shell model. Final state interactions are included within the 
distorted wave impulse approximation, by means of a relativistic optical potential, with and without 
imaginary part, and of a relativistic mean field potential. For comparison with experimental data 
the inclusive charged-current quasielastic cross section for v^- 12 C scattering in the kinematical 
conditions of the LSND experiment at Los Alamos is also presented and briefly discussed. 
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In the past few years the observation of neutrino os- 
cillations at Super-Kamiokande pj and the subsequent 
proposal and realization of new experiments, aimed at 
determining neutrino properties with high accuracy 0, 
have renovated the interest towards neutrino scattering 
on complex nuclei. In fact, neutrino detectors usually 
contain carbon or oxygen nuclei, and for a proper inter- 
pretation of the experimental results the description of 
the ^-nucleus interaction must be accurate Q . 

At intermediate neutrino energies, ranging from some 
hundreds MeV to a few GeV, ^-nucleus quasielastic scat- 
tering has been studied within several approaches Q. 
Relativistic and non relativistic studies of Random Phase 
Approximation have shown nuclear structure effects to be 
relevant only at low momentum transfers, but indications 
have been found that for future and precise data analyses 
of, e.g., atmospheric neutrino measurements, more accu- 
rate theoretical estimates may be needed. Additionally, 
very recently new attention has been drawn towards final 
state interaction (FSI) effects, which, contrary to what is 
often assumed, may still be relevant even at the relatively 
high energy = 1 GeV @. 

In this contribution we study charged-current (CC) 
neutrino-nucleus quasielastic scattering within the 
framework of a Relativistic Shell Model (RSM), already 
successfully employed to study exclusive electron scatter- 
ing and Neutral Current neutrino scattering 0. We 
compute inclusive v^- w O quasielastic cross sections for 
three values of the incident neutrino energy, namely 200 
MeV, 500 MeV, and 1 GeV, which are representative of 
the kinematical range where quasielastic scattering gives 
the main contribution to the inclusive j/-nucleus process. 

We describe the CC quasielastic scattering of neutrinos 
on a nuclear target within the Impulse Approximation 
(IA), assuming that the incident neutrino exchanges one 
vector boson with only one nucleon, which is then emit- 
ted, while the remaining (A-l) nucleons in the target are 
spectators. The nuclear current is assumed to be the sum 
of single nucleon currents, for which we employ the usual 



free nucleon expression (see yy) with the axial form factor 
parameterized as a dipole with cut-off mass Ma = 1.026 
GeV Q , and the states of the target and residual nuclei 
to be adequately described by an independent particle 
model wave function. 

To describe the bound nucleon states we use Rela- 
tivistic Shell Model wave functions, obtained as the self- 
consistent (Hartree) solutions of a Dirac equation, de- 
rived, within a Relativistic Mean Field approach, from 
a Lagrangian containing a, u> and p mesons As the 
single-particle binding energies determine the threshold 
of the cross section for every shell, in the numerical cal- 
culations, we have used the experimental values corre- 
sponding to the binding energies of the different shells. 

For the outgoing nucleon the simplest choice is to 
use plane wave spinors, i.e., no interaction is consid- 
ered between the ejected nucleon and the residual nu- 
cleus (PWIA). For a more realistic description, FSI ef- 
fects should be taken into account. In our formalism this 
is done by using distorted waves which are given as solu- 
tions of a Dirac equation containing a phenomenological 
relativistic optical potential (ROP), consisting of a real 
part, which describes the rescattering of the ejected nu- 
cleon and of an imaginary part, that accounts for the 
absorption of it into unobserved channels. In this work 
we use the ROP corresponding to the EDAD-1 single- 
nucleon parameterization presented in ref. |lfj. The use 
of this phenomenological ROP leads to an excellent agree- 
ment between theoretical calculations and data for exclu- 
sive (e, e'N) observables @; however some caution should 
be taken in extending the conclusions drawn from the 
analysis of exclusive reactions to inclusive ones. In the 
latter, unless a selection of the single nucleon knockout 
contribution is experimentally feasible, all final channels 
are included and thus the imaginary term in the optical 
potential leads to an overestimation of FSI effects. This 
is clearly illustrated in fig. 1 where we compare our theo- 
retical results with the experimental cross section corre- 
sponding to inclusive quasielastic electron scattering on 



2 




100 120 140 160 180 200 

co (MeV) 



FIG. 1: Differential cross section (da / dQ,du}) versus the en- 
ergy transfer for inclusive quasielastic electron scattering on 
12 C for a momentum transfer q ~ 400 MeV/c. The solid curve 
corresponds to the RSM with FSI described by the RMF, 
while the long-dashed curve corresponds to the real ROP and 
the short-dashed curve to the complex ROP. Finally, the dot- 
ted curve does not include FSI. Experimental data arc from 
ref. El. 



12 C for a momentum transfer equal to 400 MeV/c (sim- 
ilar results are obtained for 300 and 500 MeV/c). In our 
calculation besides the complex ROP, we also consider 
the potential obtained by setting the imaginary part of 
the ROP to zero. Additionally one may also use dis- 
torted waves which are obtained as the solutions in the 
continuum of the same Dirac equation used to describe 
the initial bound nucleon. We refer to this approach, that 
should be adequate to describe FSI at moderate energy 
transfer, as relativistic mean field (RMF). As shown in 
fig. 1, the complex ROP results clearly underestimate the 
data, while the reverse occurs for the PWIA calculation. 
On the other hand, the RMF and purely real ROP agree 
much better with experiment, particularly RMF for small 
transfer energy and real ROP for higher u. We believe 
that the RMF and real ROP results indicate a reasonable 
"band" where FSI effects should lie. 

Let us now present our results for neutrino scatter- 
ing on le O. To better illustrate our model, we start by 
neglecting FSI and comparing RSM-PWIA results with 
inclusive cross sections obtained within the relativistic 
Fermi gas (RFG) 0. This is is done in fig. [21 which 
shows the differential cross section (da/dT^) as a func- 
tion of the outgoing muon kinetic energy. With respect 
to the RFG curve calculated with no binding energy we 
observe that the RSM cross section is reduced and shifted 
towards lower values, in a way which is similar to the 
effect of an average binding energy in the RFG. In ad- 
dition the RSM cross section has a different shape, due 
to the different momentum distributions of the single nu- 
cleon shells contributing to the process. Since the various 
shells have different binding energies, the corresponding 
contributions to the cross section go to zero at differ- 
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FIG. 2: Differential cross section (da/dT^) versus the out- 
going muon kinetic energy, for the quasielastic scattering of 
muon neutrinos on 16 O and for three choices of the incident 
neutrino energy: E v — 200 MeV (upper panel), 500 MeV 
(middle) and 1 GeV (lower panel). The solid curves corre- 
spond to the RSM with no final state interaction, while the re- 
maining curves are calculated within the RFG, with fcf = 225 
MeV and e B = (dot-dashed) and e B = 20 MeV (dotted). 



ent values of T M and this gives rise to the structure of 
(da/dTfj,) observed at large T M . Fig. [21 shows that nu- 
clear model effects on the cross sections can be rather 
large at low neutrino energy, but become less relevant as 
E v increases, practically disappearing at E v — 1 GeV. 

On the other hand, the behavior of FSI effects is quite 
different, as illustrated in fig. 02 Here the results ob- 
tained with the RSM in PWIA are compared with the 
cases where FSI are described within the RMF, the real 
ROP and the complex ROP approaches. The use of real 
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FIG. 3: Same as flg|U but including FSI effects. All curves 
are calculated within the RSM model in PWIA (solid), and 
within the RMF (dashed), real ROP (dotted) and complex 
ROP (dot-dashed) approaches. 



potentials (RMF and real ROP) for describing the final 
nucleon states leads to the resonant structure observed 
for relatively high T M (that is, small energy transfer uj). 
Note that in this work we only include single-particle ex- 
citations within a mean field picture. Including residual 
interactions would make the width and number of reso- 
nances to be considerably larger. 

We observe that FSI effects produce a reduction of the 
cross section, particularly important in the case of the 
complex ROP model due to the absorption introduced by 
the imaginary term: about 60% for E v = 200 MeV and 
50% for E v = 500 MeV and E v = 1 GeV in the region 
close to the maximum. For the RMF and real ROP, the 
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FIG. 4: Integrated cross section a(E^) for the quasielastic 
scattering of muon neutrinos on 16 O as a function of the 
incident neutrino energy. The curves are calculated within 
the RFG model with ItF = 225 MeV and binding energy 
e B = (solid line) and e B = 20 MeV (dashed). The 
points correspond to RSM calculations without FSI (stars) 
and with FSI effects taken into account within the RMF 
(empty squares), real ROP (full squares) and complex ROP 
(circles) approaches. 



reduction, similar in both cases, is about 30 ~ 40% for 
E v = 200 MeV and 20% for the other energy values. 

Nuclear model effects on integrated cross sections are 
studied in fig. where the cross section cr(E v ) is plot- 
ted as a function of the incident neutrino energy. Here 
the contributions coming from the RMF and ROP res- 
onances have been included in the calculation, in order 
to respect the completeness of the set of final states pre- 
dicted by the model. These contributions are important 
at E v = 200 MeV, where they amount to about 10% of 
the integrated cross section, while at higher energies these 
effects are about 2% (500 MeV) and 1% (1 GeV). Again 
we see that within the PWIA the discrepancy between 
different nuclear models is relatively small and decreases 
with increasing neutrino energy. On the contrary FSI 
effects remain sizeable even at large E v . As in the previ- 
ous figure, the imaginary term in the ROP leads to a too 
large reduction (~ 50%) of the integrated cross section. 
The results for the RMF and real ROP models, which are 
more reliable, show a smaller, but still sizeable (~ 15% 
at E v — 1 GeV) reduction. 

The results shown in figs. Inland 0] lead us to conclude 
that FSI effects should be carefully considered in neutrino 
experiments which use oxygen based detectors. Anal- 
ogous conclusions can be drawn for the case of a 12 C 
target. 

For the purpose of comparison with the experiment, 
let us consider the inclusive 12 C(t^,/i _ )X cross section 
measured by the LSND collaboration at Los Alamos, us- 
ing a pion-Decay-in-Flight beam, with energies rang- 
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FIG. 5: Observed distribution of muon kinetic energies T M 
compared with the flux-averaged predictions of our RSM, 
in P WIA (dotted line) and including FSI within the RMF 
(solid) and purely real ROP (dashed) frameworks. The the- 
oretical distributions have been normalized to give the same 
integrated values as the experimental points, and have been 
folded in energy with a bin size of 5 MeV, the same employed 
for the experimental data. Data are from Albert et al. [Idl- 



ing from muon threshold to 300 MeV, and a large liquid 
scintillator detector ■ This experiment has been com- 
pared to very different theoretical approaches |13| ] , but an 
important discrepancy still remains. Although at these 
low energies processes different from the quasielastic nu- 
cleon knockout affect the inclusive cross section, we con- 



sider our FSI approach to be useful to estimate integrated 
cross sections. 

In fig. [5] we show the observed and calculated dis- 
tribution of events, averaged over the 1994 Los Alamos 
neutrino spectrum <p(E^) within the energy range E v = 
123.1 -j- 300 MeV. The shape and position of the maxi- 
mum of the experimental distribution are approximately 
reproduced by the three calculations, but the results that 
include FSI with the RMF potential are clearly favored 
by the data. This is consistent with the fact that includ- 
ing FSI with the mean field potential should be adequate 
at moderate kinetic energy of the ejected nucleon. How- 
ever the values we obtain for the the flux-averaged inte- 
grated cross section overestimate the measured cross sec- 
tion by approximately 50%. More precisely, in the RSM, 
we obtain (a) = 20.5 (PWIA), 16.8 (RMF) and 15.1 (real 
ROP) 10~ 40 cm 2 . Additional corrections due to the out- 
going muon Coulomb distortion, evaluated within the ef- 
fective momentum approximation |14|. further increase 
these numbers by 5 4- 10%. The corresponding final mea- 
sured experimental value is 10.6±0.3±1.8 10~ 40 cm 2 [l2| . 
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